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Introduction 

Until  the  late  80's  information  on  Soviet  underground  nuclear 
explosions,  let  alone  seismic  recording  on  Soviet  soil,  were  not  generally 
available  to  the  seismological  community.  The  dramatic  changes  since  then 
have  resulted  in  access  to  a  wide  variety  of  seismic  data  from  Soviet 
explosions.  In  this  report  we  present  a  data  set,  unique  to  the  field  of 
seismic  verification.  The  data  includes  a  large  number  of  waveforms  from 
in-country  seismological  station  Borovoye  for  the  Soviet  peaceful  nuclear 
explosions  with  announced  yields  and  origin  times,  some  with  physical 
and  mechanical  conditions  at  the  test  site.  The  waveforms  were  recorded 
by  digital  station  of  different  types.  In  this  report  we  have  summarized  and 
rewiewed  information  -on  122  explosions,  and  55  waveforms  in  this  data 
set,  which  contains  recordings  obtained  in  the  course  of  more  than  15 
years.  As  the  characteristics  of  recording  instruments  were  changing 
during  this  period  of  time  available  information  on  instrument^ 
calibrations  is  also  described  and  reviewed. 

We  also  have  attempted  to  describe  the  observed  pecularities  of  P- 
wave  by  using  a  simple  source  function  and  magnitude  correction 
which  take  in  attention  such  parameters  as  the  velocity  of  longitudinal 
wave,  density,  moisture,  gase  content  and  depth  of  explosion.  More  detail 
analysis  was  conducted  for  the  explosions  conducted  in  salt  as  there  is 
representative  sample  and  they  were  observed  on  local  distances. 


Peaceful  nuclear  explosions 

122  explosions  were  conducted  in  USSR  in  the  period  from  1965  to 
1988.  The  map  in  Figure  1  shows  the  location  of  explosions  sites,  which 
are  squares  and  seismological  station  Borovoye  which  is  asterisk.  The  size 
of  square  scaled  by  magnitude  of  explosion.  The  boundaryes  of  major 
tectonic  elements  of  the  USSR  are  portrayed  with  lines  in  Figure  1  showing 
that  the  parts  traverse  a  diversity  of  tectonic  structures. 

Information  on  date,  origin  time,  location  and  body  wave 
magnitude  m.,  reported  by  ISC,  is  compiled  in  Table  1.  Precise  times, 
based  on  nonseismological  information,  have  been  collected  for  78 
explosions  (Nuclear  Explosions  in  the  USSR,  1994)  .  The  epicenter 
estimates,  entered  in  Table  1,  are  limited  to  those  of  ISC  as  they  are  the  only 
type  determination  consistently  applied  all  the  explosions,  with  exception  for 
15  -  explosions  which  location  based  on  nonseismological  information  . 
(marked  by  *)  have  published  by  Sultanov  et  all  (1993).  To  estimate  the 
errors  in  locations  and  origin  times  ,  obtained  at  ISC,  we  have  compared 
nonseismological  and  seismological  information.  We  have  recieved  the 
systematic  bias  about  2.6±0.4  sec  which  agrees  with  the  estimation  of  the 
arrival  time  curve  of  P-wave  by  (Kogan  S.D.,1976).  The  root  mean  square 
error  in  location  has  been  estimated  as  6±4  km.  As  results  the  origin  tunes 
for  44  explosions  ,  based  on  seismological  data  ,  are  the  ISC  origin  times 
corrected  for  the  systematic  bias. 

Table  2  summurizes  the  information  on  the  type  of  environmental 
media  contained  explosions  and  summerizes  the  quantitive  value  such 
parameters  of  media  as  the  longitudinal  wave,v  ,  density, p,  gaseous 
content, y  ,and  moisture,®.  The  information  on  yield, q,  and  reduced  depth 
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is  also  compiled  in  Table  2.  The  velocity  ,  entered  in  Table  2  ,  was 
estimated  by  seismic  and  acustic  sounding.  The  density  was  estimated  from 
the  sample  of  rock.  The  relative  moisture  is  the  part  of  water  which  was 
vaporized  from  the  sample  of  rock  ,  heated  up  to  100°  C.  The  relative 
gaseous  content  is  the  number  of  gases  which  was  recieved  from  the  sample 
of  rock  ,  heated  up  to  1000°  C.  In  that  cases  when  the  environmental  media 
was  unhomogenity  we  used  the  parameter  as  the  average  for  the  layer 
which  thick  is  equal  two  caveme  radius  rc.  That  is 

a  =  £  ht  and  £  /»,  =  2rc  where 

ai  h  i  are  the  value  parameter  of  i  in  the  layer  with  the  thick  of  h  i.  The  table 
shows  that  the  parameters  are  limited  1.8  <  v  ^4.4  km/sec,  170  ^h 
<  2500  m,  1.7^  p<  3.0  g/cm3,  0<a><30%,  0<y<40%,  300<A<3300  km. 

Of  all  the  explosions  with  known  conditions  at  the  test  site  until  we  have 
collected: 

53  records  with  P-wave, 

25  records  with  Lg  wave, 

7  records  with  surface  wave. 

Plot  record  section  of  vertical  component  for  some  explosions  is  shown 
on  the  Figure  2.  The  analysis  of  these  records  shows  that  from  the  distances 
about  1960  km  the  maximum  amplitude  in  the  S-  and  Lg-  waves  becomes 
much  lower  than  the  P-wave  amplitude.  In  the  range  of  epicentral  distances 
from  300  km  to  1600  km  the  maximum  amplitude  of  S-  and  Lg-waves  is  a 
comparable  with  P-wave  amplitude.  The  exception  is  the  set  of  explosions  in 
the  Uralsk  region,  being  carryed  out  in  salt,  at  the  distance  about  1160  km. 
The  amplitudes  of  S-  and  Lg-waves  were  lower  than  for  P-wave  ones  in 
about  3-4  times.  It  is  worth  mentioning  that  nonseismological  information  on 
the  yield  or  the  conditions  at  the  test  site  of  explosion,  carryed  out  on 
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2.4/09/83  (origin  time  is  05h09m57.4  sec),  shown  in  Table  1  are  not  seems  to 
be  correct.  In  the  same  conditions  of  media  where  this  set  was  carryed  and 
similarity  of  wavefonn  of  all  six  explosions,  the  maximum  amplitude  of  P- 
wave  differs  a  great  deal  from  the  amplitude  for  other  explosions.  The 
magnitude,  derived  from  the  data  of  in-country  stations  ,  for  the  explosion 
under  question  is  4.85  while  for  other  5  explosions  is  in  the  range  5.18-5.19. 
It  follows  the  the  yield  is  3.9  kt  according  to  formula  (12)  below.  The  relative 
yield  estimation  on  Borovoye  data ,  presented  in  (Sultanov  D.J.  et  all,  1993  ) 
is  3.4.  That  is  the  yield  is  about  2.5  kt. 


The  Borovoye  Geophysical  Observatory 

The  Borovoye  Geophysical  Observatory  is  one  of  the  oldest 
seismic  stations  of  the  world  in  terms  of  the  digital  seismic  recording.  The 
observatory  has  the  richest  archive,  comprizing  the  digital  material  of 
various  lands  of  recording  systems  working  from  1965  until  nowadays. 
There  is  therefore  the  Borovoye  station  draw  the  attention  of  the  world 
community  of  scientists. 


Station  characteristics 

The  Borovoye  seismic  station  (BRVK)  with  coordinates  53°.083  North, 
70°.250  East  is  located  in  the  North  Kazakstan  about  70  km  to  the  South 
from  Kokchetav  town. 

In  terms  of  geology  the  Borovoye  region  presents  as  a  flat  denudation 
country  and  hills  with  altitude  change  from  160  to  500  meters.  The  region 
is  made  up  of  Paleozoic  and  prior-Paleozoic  dislocated  rocks, 
breaked  by  intrusive  igneous  rocks  of  various  age.  On  the  surface  there  are 
loose  stuff,  presented  by  loosened  granite,  sands,  loamy  soils,  clays  the 
thick  of  which  may  vary  from  severel  meters  to  several  hundreds  meters.  The 
P-wave  velocity  in  the  layers  varies  from  0.65  to  2.8  km/sec.  From  seismic 
studies  it  follows  that  the  thick  of  granite  layer  is  estimated  as  8-10  km 
and  basalt  layer  -  as  35  km.  The  thick  of  crust  vary  from  40  to  50  km. 

The  main  tectonic  element  in  the  block  structure  of  the  tectonic 
framework  of  the  region  is  extremely  complicated  by  the  faults  of  the 
first  order,  they  having  north-west  direction.  The  faults  of  the  second 
order  have  different  directions.  From  gravity  and  magnetic  studies  it 
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follows  that  the  faults  of  the  first  order  penetrate  through  the  granite  layer  as  a 
whole. 

At  the  station  place  the  granite  layer  comes  out  on  the  surface  and 
almost  monolithic  granite  are  on  the  deep  to  vary  from  6  to  10  meters. 

The  tectonic  structure  of  the  North  Kazakhstan  is  known  from  the  deep 
seismic  studies  which  was  earned  out  by  means  of  numerous  profiles 
(Antonenko  A.,N.,1984).  The  main  profile  was  the  profile  of  extent  about 
1000  km  from  Temirtau  to  Petrorpavlovsk  (Figure  3). 

Side  by  side  with  deep  faults  the  massifes  of  rocks  of  Agnostozoic  era 
are  the  main  tectonic  element  of  North  Kazakhstan  structure.  The 
Borovoye  station  is  situated  on  one  of  these  massifes,  named  Kokchetav 
anticlinorium.  Note  should  be  taken  that  the  crust  was  studied  in  detail  in  the 
profile  zones.  The  profile  observation  system  of  seismic  waves  have  given 
the  possibility  to  record  the  waves  refracted  on  Moho  and  on  the 
discontinuties  in  the  crust.  Thus  the  refracted  waves  from  Moho 
discontinuty  were  well  observed  along  the  profiles  on  large  distances. 
However,  beneath  the  region  of  Borovoye  the  sharp  eliminate  of  amplitudes 
of  P-waves  refracted  on  Moho  discontinuty  was  being  observed. 
Dissapearance  of  the  waves  refracted  on  Moho  may  be  explained  by 
existance  of  strongly  dislocated  zones  in  the  mantle  beneath  Borovoye 
region. 

The  second  main  rezult,  carried  out  by  means  of  deep  seismic  study 
was  that  beneath  Borovoye  region  the  crust  discontinuties  were  not 
revealed.  So  the  crust  under  the  place  may  be  considered  as  "transparent" 

Until  1972  the  seismometers  of  the  station  were  installed  in  the  gallery  in 
weekly  cracked  granite  with  absolute  altitude  of  about  340  meters  and 
from  1972  they  were  installed  in  the  mine  with  deep  of  15  meters  from  the 
surface  and  altitude  to  be  equal  315  meters  in  practically  monolithe  granite. 
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The  seismic  station  was  set  up  for  the  observation  of  seismic  waves 
generated  by  nuclear  explosions  at  the  Nevada  test  site.  The  first  records  of 
seismic  waves  in  Borovoye  showed  that  the  signals  from  Nevada  were  two 
or  three  times  greater  than  would  be  expected. 

While  developing  of  seismic  recording  several  kinds  of  system  of 
registration  were  used.  The  first  digital  system,  known  as  KOD,  began 
recording  in  1966  and  continuous  worked  from  1967  to  1973 
(AdushkinV.,V.,and  V.,A.  An,  1992).  It  was  based  on  three  component 
short-period  seismometers  and  was  one  of  the  oldest  digital  system,  operating 
in  the  world  in  1960s  and  1970s.  Later,  other  digital  systems  were 
installed.  They  are  known  as  STR-SS  and  STR-TSG.  The  digital  system 
STR-SS  was  intended  mainly  for  low-gain  recording  and  the  system  STR- 
TSG,  which  operated  with  twenty  four  channels,  most  recorded  a  two  gain 
levels.  STR-TSG  system  is  based  on  long-period  and  short  period  Kimos 
seismometers.  Both  STR-SS  and  STR-TSG  have  been  operating  from  1973 
up  to  the  present.  The  dynamic  range  of  every  channel  makes  up  60 
dB.  The  Figure  4  demonstrates  amplitude-frequency  response  of  some 
digital  channels,  operating  at  Borovoye  station. 

The  seismic  station  is  situated  at  the  place  of  low  level  of  noise  to  be 
both  natural  and  industrial.  The  amplitude  of  the  noise  of  short  period  range 
vary  from  one  to  ten  nanometers  and  the  noise  of  long  period  range 
composes  r00-200  nm. 


MORFOLOGY  of  SEISMIC  SIGNALS 


General  description 

When  P-waves  propagate  through  the  Earth  the  high  frequencies  are 
preferetially  attenuated  so  that  individual  pulse  are  smeared  out  and  may 
overlap.  Also,  as  the  recording  system  has  limited  pass  band,  signals  are 
further  distorted  on  recording.  Analysis  and  interpretation  of  short-period 
seismograms  may  thus  be  difficult.  The  regional  observed  waveform  is  very 
complicate  due  to  the  path  of  propagation.  There  are  a  number  of  seismic 
phases  (reflected,  refracted,  converted),  generated  on  the  stucture 
inhomogenities  in  the  crust  and  upper  mantle.  This  is  a  good  aid  to 
interpreting  short-period  explosion  recordings  is  to  compute  seismograms 
using  model  of  the  Earth,  source  and  seismograph  and  try  and  match  these 
the  observed. 

Comparing  waveform  of  P-wave  from  PNE  at  the  distances  300- 
3300  km  reveals  their  substantial  discrepency.  Nevetheless,  it  takes  place 
an  extremal  similarity  that  could  be  presented  as  a  similar  initial  part  of  P- 
waves,  which  is  consisted  of  three  phases  that  are  clear.  At  the  distances 
up  to  1000  km  the  first  arrival  has  a  period  of  0.5-0. 6  sec,  the  period  of  the 
second  phase  is  1.0-1. 2  sec,  then  again  the  phase  with  period  about  0.5- 
0.6  sec  appears.  It  is  that  the  maximum  amplitude  is  associated.  At  the 
distances  about  1700  km  another  picture  can  be  observed:  the  first  arrival 
has  a  larger  period,  the  second  has  a  smaller  one  and  the  third  phase  has  a 
larger  period  again.  The  example  of  record  is  shown  on  the  Figure  4.  This 
pecularity  is  typical  for  events  on  the  regional  distances. 

At  the  distances  up  to  1000  km  the  first  arrivals  is  a  wave  with  the 
velocity  about  8.15  km/sec  associated  with  Moho  discontinuity.  At  the  far 
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distances  the  first  arrival  seems  a  head  wave  propagating  in  the  upper  mantle 
with  the  velocity  about  8.48  km/sec  and  associated  with  discontinuity  on  the 
depth  about  150-200  km. 

However,  there  is  the  substantial  pecularity  of  records  of  the  events 
carried  out  to  the  North  of  the  station.  The  maximum  amplitude  in  P-wave  is 
at  the  later  time  than  for  the  signals  arriving  from  other  directions.  The 
maximum  amplitude  arrival  time,  Pmtime,  and  the  first  arrival  time,  Ptime, 
summerizes  in  the  Table  3  and  the  examples  of  seismogram  shown  on  the 
Figure  6  (a,b,c).  It  can  see  that  the  time  delay  of  the  maximum  amplitude 
arrival  relative  the  first  time  arrival  is  about  40  sec  (Figure  6  a,b).  The 
parameters  of  explosions  ,  shown  on  the  Figure  6,  presented  in  Table  4.  One 
can  assume,  that  observed  pecularity  could  be  related  with  the  disturbance  of 
Moho  discontinuity  to  the  North  from  Borovoye  station,  which  was  detected 
by  deep  seismic  sounding,  or  with  a  vast  fault  which,  as  supposed,  penetrates 
to  the  upper  mantle  (Antonenko,  1972). 


Incident  angle  and  azimuth 

Because  of  strong  interference  in  the  wave-train  in  the  5  sec  interval 
after  the  first  arrival ,  which  contains,  possibly,  two  arrivals  on  the  distance  of 
1200  km  and  three  arrivals  on  the  distance  1700  km,  detailed  study  was 
carryed  out  of  these  arrivals. 

First  of  all  the  incident  angle  and  azimuth  of  coming  waves  were  estimated  by 
using  the  analysis  of  eigenvalue  and  eigenvector  of  the  covariance  matrix  in 
the  moving  time  window.  For  instantanous  motion  ut  (i=l,2,3),  the  cross 
power  matrix  is 


(1) 


1  N 

C*/  =  T7SW/  k,l=l,2,3 
**  1  =  1 

The  eigenvalue  X  and  eigenvector  v  of  matrix  C  can  be  found  from 
(C  -  X2I)vk  =  0,  m=l,2,3  (2) 

and  are  ordered  by  size,  I  is  the  unit  matrix.  The  normalized  eigenvalue 
v=(vlm ,  v2ra  ,v3m )  corresponding  to  the  largest  eigenvalue  X  may  than  be  used 
as  an  estimate  of  the  dominant  signal  direction.  The  apperent  azimuth  and 
angle  of  incidence  are  to  be  found  from 

tan  Az  =  — ,  Cost  =  vl2 .  (3) 

Vi„ 

The  size  of  X  can  be  used  as  a  measure  of  the  stability  of  this  procedure  since 
it  will  be  a  unit  for  purely  rectilinear  motion.  These  directional  parameters  are 
calculated  from  the  particle  motion  of  the  code,  hence  they  do  not  necessarily 
give  the  actual  direction  of  propagation.  The  polarization  parameters  as 
rectilinearity  was  also  calculated.  The  measures  of  rectilinearity  are  based  on 
the  eigenvalue  of  covariance  matrix  employed  in  the  power  technique: 

red  =  1  -  (A3  +  X2)  /  2Aj,  (4) 

where  X\  ,Xi  ,Xi  are  the  eigenvalues  (largest  to  smallest)  of  cross  power 
matrix  C.  hi  Figure  7  we  display  the  wave  field  measured  introduced  above. 
For  the  P-waves,  the  estimated  azimuth  is  close  to  the  source-receiver 
azimuth  of  274°  for  the  explosions  on  the  distance  about  1160  km  and  is 
close  to  the  source-receiver  azimuth  of  255°  for  the  explosions  on  the 
distance  of  1715  km.  The  azimuth  is  stable  for  the  interval  of  about  6  sec 
relatively  the  first  arrival.  The  azimuth ,  incident  angle  and  the  ratio  of  middle 
Xi  and  maximun  eigenvalue  X\  related  to  different  phase  arrival  are  presented 
in  Tables  5,6.  One  can  see  from  Tables  5,6  that  there  are  common  features  in 
behavious  of  measure  of  wavefield  parameters.  The  apparent  azimuth  is  fairly 
stable  for  the  first  6  sec  of  the  seismogram,  but  some  fluctuations  are  seen 
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and  the  systematic  azimuth  deviation  with  increasing  time  too.  Effect  of 
decreasing  azimuth  with  time  is  observed.  The  apparent  angle  of  incidence 
shows  considerable  variability  but  clearly  indicates  on  the  distance  of  1716 
km  on  the  increasing  in  the  begining  of  signal  and  the  decreasing  in  the  end 
of  signal.  The  maximun  value  is  44.5°  .  On  the  distance  of  1160  km  the 
incidence  angle  decreases  with  time. 

High  rectiliniarity  values  are  generally  associated  with  stable  azimuth 
values,  and  at  the  places  where  the  azimuth  fluctuates,  the  values  indicate  the 
waves  are  not  rectilinear.  We  have  done  this  analysis  in  order  to  try  to  define 
the  wave  types  which  can  appear  on  the  seismogram. 

Summarizing  the  results  of  above  analysis  it  is  possible  to  conclude  that  there 
are  several  linear-polarized  signals  in  the  first  arrival  of  P-wave  with  different 
velocities.  It  is  usual  for  the  distance  of  1716  km.  We  can  suggest  that  these 
pecularities  are  related  to  the  reflected  and  refracted  waves  in  the  upper 
mantle  .  hi  particular  the  first  arrival  on  the  distance  up  to  1000  kilometers  is  • 
the  head  wave  Pn  due  to  Moho  discontinuty,  and  the  following  phase  due  to  a 
more  deep  part  of  the  upper  mantle.  On  the  distance  up  to  1000  kilometers 
the  Pn  phase  marked  the  arrival  of  longitudinal  waves  traveling  in  the 
uppermost  mantle.  It  is  known  from  numerous  investigations  that  the  travel 
time  curve  of  Pn  is  practically  a  straight  line.  The  slope  has  been  found  about 
13.5  sec/degree,  which  corresponds  to  wave  velocity  about  8.2  km/sec 
(Passechnik ,  1972)  .  The  first  arrival  on  the  distance  1716  km  do  not  mark  as 
P  wave  refracted  in  gradient  velocity  of  the  earth  because  of  the  complex 
waveform.  We  suppose  that  the  first  arrival  can  be  a  head  wave  generated  by 
the  discontinuity  on  the  depth  150-200  km.  It  can  be  supposed  that  there  is 
low  velocity  layer  on  this  depth  because  the  great  amplitude  on  the  2-3  sec 
after  the  first  arrival.  The  maximum  amplitude  arrival  seems  to  be  reflected 
wave  from  discontinuty  on  the  depth  of  150-200  km. 
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On  the  distance  above  2100  km  the  waveform  is  very  simple  and  we 
can  suppose  that  the  P-wave  is  the  first  arrival  because  of  the  longer  period 
and  faster  velocity. 
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Source  model  and  scaling  low 


Iii  the  vibrating  sphere  problem,  an  explosion  is  modelled  at  and 
beyond  some  critical  distance  -  elastic  radius,  where  the  material  behaves 
elastically  by  a  radial  stress  applied  uniformaly  over  a  spherical  surface  . 
The  spherical  surface  which  separates  inelastic  from  elastic  response 
has  been  called  by  equivalent  cavity  by  Sharpe  (  1942). 

For  a  step  in  pressure  p(t)=PoH(t) , where  P0  is  the  amplitude  of  the  pressure 
applied  to  the  inside  of  a  spherical  surface  at  r=Ro  and  H(t)  is  Heaviside 
step  function  that  is  H(t)=l  for  t>0  and  H(t)=0  for  t<0. 

In  this  case 


where  vp  and  vs  are  longitudinal  and  shear  velocity  ,  n  is  shear  modulus, 
and  (o  and  kp  are  angular  frequency  and  wavenumber ,  respectively. 
Spectra  of  the  direct  P-wave  have  the  obvious  asymtotic  behaviour 

—  for  (6)  and 

*  4  vpfi  r 


«p(r,e>)  = 

P  4  vp/t  R0co 


for  o)  >  o)o 
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With  the  base  on  these  relations  one  can  get  the  comer  frequency  co  o 
as  a  point  of  crossing  of  the  asymtotes  that  is  the  point  where  high-frequency 
asymptote  crosses  low-frequency  one. 

To  use  these  relations  for  the  explosions  conducted  in  the  media  with 
different  mechanical  and  physical  behaviour  we  need  the  knowledge  what 
way  PqRI  depends  on  above  parameters.  The  work  (Laushkin  V.A.  and 
Ovtchinnikov  V.M.,  1992)  offers  a  simple  model  for  taking  into  account  a 
physical  and  mechanical  features  of  media  when  estimating  the  yield  of 
explosion.  It  was  done  by  introducing  the  magnitude  correction  according  to 
the  following  relation 

Am  =  1.5 / 1.75  log  —  +  0.39  log—  +  31og(l  +  2.9 <y)  +  31og(l  +  0.62^)  +  A mh  ]  (8) 

\  Pa  ' 


A mh  =  0.32 log  =  for  hlh<  5 

(9) 

and  Amh  =  0.321og=-+0.4  for  h/h<  5, 

(10) 

where  h  =  100 \fq  . 

In  this  relation  v  is  the  longitudinal  velocity,  p  is  density  of  media 
contained  explosion  ,  to  is  relative  moisture  ,  y  is  relative  gase  content,  h  is  a 
depth  of  and  q  is  yield  of  explosion.  Index  0  refers  to  some  reference  media 
for  which  the  relation  between  m  and  q  is  known.  Figure  8  presents  the 
relation  m=f(q)  for  the  explosions  from  Table  1  and  magnitude  correcton  from 
Table  2.  On  this  picture  (a)  we  can  see  a  substantial  scatter  of  experimental 
data  relatively  the  line 

m=0.771ogq  +4.4.  (11) 


( 


Introducing  die  magnitude  correction  ,  calculated  in  according  to  above 
relation  ,  we  can  see  from  picture  (b)  that  the  experimental  points  place  more 
close  relatively  the  line 

m-Arn=0.811ogq+4.4.  (12) 

The  relative  root  mean  square  error  is  20%. 

To  give  a  physical  interpretation  of  magnitude  correction  A m  we  will 
use  the  Sharpe  model  as  die  the  source  function. 

Following  Nersesov  et  all  ( 1991)  let  the  changing  of  elastic  energy  of 
media  be  proportional  to  the  yield  of  explosion,  e.g.  AU=Cq,  where  C 
depends  upon  the  characteristics  of  media  ,  containing  the  explosion.  On  the 
other  hand,  the  work  made  on  the  elastic  radius  by  pressure  Po  is 


f  ar\ 

AU  =  AW  =  f  4r2P0ndr  =  4P0R%  —hr,  (13) 

i  W 

where  AR  is  motion  of  elastic  boundary  due  to  pressure  P0  on  R0.  From  the 
condition  that  the  pressure  on  R0  is  equal  to  the  elastic  limit  of  media  ,Y  we 
can  write 


(14) 


where  ks  is  adiabatic  volumetric  modulus. 

Therefore  ,  from  relations  (13),  (14)  we  receive 

PX=Kcq  ,  (15) 

where  K  =^~  is  constant  of  nonlinear  coupling  of.  explosion,  which 
c  4  Yn 

depends  on  the  elastic  parameters  of  media  ,  porosity,  moisture  and  gase 


content. 


For  the  spectral  amplitude  for  the  regional  (teleseismic  )  distances  by 
using  the  Sharpe  model  as  the  spectral  source  function  we  can  write 


U(A,a>)  =  ^qG(A)Qxp(i7rftt) 


1- 


f  V 
CO  ' 


+  /- 


.2(0 

0>n 


(16) 


where  Ke  is  constant  of  elastic  coupling,  G( A)  is  the  geometries  spreading 
fimetion,  and  t*  is  the  apparent  attenuation. 

Body  wave  magnitude  for  the  period  of  1  sec  is 


m  =  log|w|  +  i(A) . 

To  take  in  attention  the  low-frequency  asymptote  it  follows 

f  g  'N 

m  =  log  q  +  log^-J  +  b(  A). , 


(17) 


where  b{ A)  is  the  calibrate  function  (Vanek  U,  Kondorskaiy  N.,V.,  et  all, 
1962)  that  is  the  correction  for  the  attenuation  on  the  path  of  propagation  and 
the  geometries  spreading. 


Thus  magnitude  correction,  introduced  above,  A m  =  log 


Spectra  of  some  explosions  in  salt 

The  instrument  -corrected  amplitude  spectra  can  be  parametrized  as 
U(co,A)  =  S(co)G(A)Qxp(-^-)  (18) 

where  5(0)  is  the  source  spectrum,  G( A  )  is  the  geometric  spreading 
function  and  last  term  is  the  apparent  attenuation  for  the  travel  time  t.  The 
geometric  spreading  of  head  wave  is  a  complicated  function  of  the  velocity 
gradient  in  mantle  and  it  is  probably  frequency-dependent.  Therefore,  a 
simple  parametrization  such  as  (18)  is  not  likely  to  be  applicable  to  a  head 
wave  for  a  broad  frequency  band  and  distances  range.  To  exclude  the 
influence  a  geometric  spreading  function  we  have  compared  four 
explosions  conducted  in  salt  in  close  distance  range. 
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Spectra  were  calculated  for  two  time  windows  which  duration  was  2 
sec.  The  begining  time  of  windows  relatively  the  first  arrival  time  were  0 
sec  and  3  sec  consequently. 

P-wave  spectra  of  four  explosions  with  yield  103,  58,  23,  and  9.3  kt 
in  salt  ,  corrected  for  instrument  responce  and  for  attenuation  (Q-factor  is 
equal  1200)  are  shown  in  Figures  9,10.  We  assumed  that  the  corner 
frequency  of  these  explosions  should  be  below  4  Hz,  and  the  low- 
frequency  asymptote  is  proportional  co1  for  the  signal  spectra  of  head  wave 
(Zvolinskiy,1958)  and  is  the  plate  for  the  signal  spectra  for  the  second  time 
windows.  The  high-frequency  asymptote  is  proportional  co”  where  r|=3  for 
ahead  wave  and  r|=2  for  the  second  time  window. 

We  can  see  from  Figure  9  that  the  spectra  for  the  first  time  window 
were  very  slightly  modulated  by  periodic  function.  In  the  low  frequency 
range  the  specta  fall  down  as  co1  .  Spectra  of  the  second  time  window  on 
Figure  10  have  distinct  relative  minimum,  and  is  the  plate  in  the  low- 
frequency  range.  Figures  9  and  10  show  also  the  theoretical  far-field  . 
amplitude  spectra  of  the  Sharpe  model  calculated  for  explosion  with  yield 
of  100  kt  following  Evernden  et  all  (1986).  The  corrected. spectra  of  as  the 
first  time  window  as  the  second  time  window  doesn’t  show  a  clear  corner 
frequency,  perhaps,  because  of  the  complicated  multiple  arrival.  However 
it  should  be  emphasized  that  the  low-frequency  asymptote  is  in  the 
agreement  with  our  assumptions.  But  we  can  not  explain  observed  relative 
minimun  in  spectra,  as  usual,  by  influence  of  the  depth  of  events. 
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Conclusion 


1.  There  was  collected  seismological  and  nonseismological 
information  on  122  peaceful  underground  nuclear  explosions  on  the 
territory  of  the  USSR  which  includes  the  data  on  date,  origin  time, 
coordinates,  depth,  yield,  magnitude  as  well  as  data  on  parameters  of  media, 
where  the  explosions  were  carried  out,  such  as  velocity  of  the  longitudinal 
waves,  density,  moisture,  gaseous  content. 

2.  53  digital  short-period  records  of  Borovoye  seismic  station  situated 
in  the  North  Kazakhstan  were  collected,  edited,  and  analysed.  The  channel 
responces  accompany  these  records. 

3.  The  analysis  of  prominent  pecularities  of  waveforms  allows  to 
assume  that  on  the  distances  up  to  1000  km  first  arrival  is  a  head  wave 

generated  by  Moho,  and  on  the  distances  in  the  range  1200  -1900  km  first 

,  ¥ 

arrival  is  a  head  wave  associates  with  the  discontinuity  in  upper  mantle  at 
the  depth  about  150-200  km. 

4.  The  influence  of  media  features  on  the  seismic  efficiency  was 
investigated,  and  the  relation  between  the  characteristics  of  media 
parameters  of  explosion  on  the  base  of  Sharpe  model  was  established. 

5.  The  explosions  carryed  out  to  the  North  from  the  seismic  station 
have  a  specific  waveform  due  to,  probably,  the  disturbance  of  Moho 
discontinuity. 

6.  The  spectra  of  explosions  in  salt,  corrected  for  instrument  responce 
and  for  attenuation,  don’t  show  a  clear  corner  frequency  ,  perhaps  because 
of  the  complicated  multiple  arrival. 
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Tables 


Table  1.  The  list  of  peaceful  underground  nuclear  explosions  conducted  in  the 
USSR  with  date,  type  of  media  contained  explosion,  co-ordinates, 
origin  time,  depth,  yield,  and  ISC  magnitude. 

Table  2.  The  list  of  peaceful  underground  nuclear  explosions  conducted  in  the 
USSR  with  yield ,  Hreduced  =h/h  (h  is  depth,  and  h  =100q1/3  ) 
magnitude  correction  for  the  depth,  Mcor,  physical  parameters  of 
media,  such  as  velocity,  density,  moisture  ,  gaseous  content,  magnitude 
correction  according  to  formula  (8,  page  1 5),  and  ISC  magnitude,  m. . 

Table  3.  The  results  of  measurements  of  the  travel  time  of  the  first  arrival , 
Ptime ,  and  of  the  maximum  amplitude  arrival,  Pmtime,  maximum 
amplitude  ,  Amax  in  nanometers,  and  period  of  maximum  amplitude,  T. 

Table  4.  The  basic  parameters  of  three  pairs  of  explosion. 

Table  5.  Backazimuth,  incident  angles,  and  relation  of  eigenvalues  for  the  first 
explosion ,  conducted  1983,  September,  24. 

Table  .  Backazimuth,  incident  angles,  and  relation  of  eigenvalues  for  the  first 
explosion ,  conducted  1984,  July,  21. 
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Figure  captures 


Figure  1.  Map  of  peaceful  nuclear  explosions  in  the  USSR.  Events  shown 
square  with  size  scaled  by  magnitule.  The  Borovoye  station 
marked  by  asterisk. 

Figure  2.  Plot-record  section  of  vertical  component  displacement  for  some 
explosions  in  the  distance  range  from  1600  to  2000  km.  The  time 
reduced  to  the  first  arrival  of  each  explosions.  A  row  marks  the 
seismic  phases. 

Figure  3.  The  main  profile  of  deep  seismic  sounding  in  the  North  Kazakhstan. 
The  lines  show  the  different  parts  of  profile  and  the  close  circles 
show  shot-points,  and  triangle  is  Borovoye  station.  PI,  P2  are 
presumed  head  wave  and  reflected  wave  generated  on  the 
discontinuity  in  upper  mantle  at  the  depth  about  150-200  km. 

Figure  4.  The  typical  amplitude-frequency  responce  of  digital  channels  at 

Borovoye  station,  based  on  KS,  KSM,  DS,  and  DSM  seismometers. 

Figure  5.  The  example  of  record  and  explosion  carried  out  in  salt  at  the 
distance  about  1700  km  .  The  rows  show  the  arrival  of  different 
phases. 

Figure  6.  The  comparison  of  seismograms  observed  from  explosions  carried 
out  to  the  North  from  Borovoye  station  with  seismogram  from 
explosions  from  another  direction:  distance  about  850  km  (a), 
distance  about  990  km  (b),  and  distance  about  1720  km  (c). 

Figure  7.  The  examples  of  the  polarization  characteristics  of  P-waves  from  ' 
explosions  in  salt  at  the  distance  of  1160  km  (a),  and  of  1715  km 
(b):  trace  1  is  recording  of  vertical  component  of  motion,  trace  2  is 
the  biggest  eigenvalue  trace  3  is  the  relation  of  eigenvalue  XJXz, 
trace  4  is  backazimuth,  and  trace  5  is  incident  angle. 

Figure  8.  The  dependence  of  magnitude,  m  from  yield  of  explosion,  q  without 
madnitude  correction  (a)  and  with  magnitude  correction  (b).  The 
points  are  experimental  data  from  Table  1 . 

Figure  9.  Amplitude  spectra  of  displacement  of  initial  2  sec  of  P-wave  of  four 
explosions  in  salt  with  yield:  103  kt  (a),  58  kt  (b),  23  kt  (c),  and 
9.3  kt  (d).  Theoretical  spectrum  for  yield  100  kt  according  to 
Evemden  et  all  is  presented  by  the  line  (e).  Dates  of  explosions 
placed  in  the  right  comer  of  picture. 

Figure  10.  Amplitude  spectra  of  displacement  for  the  window  from  3  to  5  sec 
relatively  the  first  arrival  time  .  Denotes  are  the  same  as  at  Figure  9. 


f  i 


40 


Figure  1 


o 


Time  (sec)  [VM  OFF] 


Figure  2 


42 


t 


Figure  3 


uojoiuu/^unoo 


46 


ie,sec 


uojoiw/}unoo 

48 


t 


time, sec 
Figure  6  c) 


zm  7  ~  ^9+01 X  mu/)[M03 


49 


f 


ime.sec 

Figure  7  a) 


time, sec 
Figure  7  b 


■7907-t4-z?,ffes3 


